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Abstract 

Although f{R) modification of late time cosmology is successful in explaining present 
cosmic acceleration, it is difficult to satisfy the fifth-force constraint simultaneously. Even 
when the fifth-force constraint is satisfied, the effective scalar degree of freedom may move 
to a point (close to its potential minima) in the field space where the Ricci scalar diverges. 
We elucidate this point further with a specific example of f{R) gravity that incorporates 
several viable f{R) gravity models in the literature. In particular, we show that the 
nonlinear evolution of the scalar field in pressureless contracting dust can easily lead to 
the curvature singularity, making this theory unviable. 
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1 Introduction 


Modern cosmological models suffer from a major theoretical difficulty, namely, the dark energy 
problem. The problem stems from the high-precision observational data with strong evidence 
that the universe is undergoing a phase of accelerated expansion in recent times [I]. Many 
models have been constructed for an explanation of the late time acceleration of the Universe, 
and those models can be mainly divided into two categories: dark energy models that change 
the matter content of the Universe and modihed gravity models that alter the Einstein gravity. 

Dark energy is generally modelled by the vacuum energy with an equation of state parameter 
w = —1 or by a dynamical scalar held (dubbed as quintessence) with tc ~ — 1 [2]. Even though 
both the cosmological constant and the quintessence held are observationally consistent in 
explaining the cosmic acceleration, the origin of these matter sources from fundamental physics 
is not well-understood [3]. 

On the other hand, late time acceleration can also be obtained from the modihcations 
of General Relativity, which introduces an extra degree of freedom to the gravitational sector 
itself. The simplest model in this category are f{R) theories, where f{R) is an arbitrary (usually 
analytic) function of R with R being the Ricci scalar laiSlEllT]. The equation of motion of 
the extra degree of freedom other than the graviton is of second order [8], and the models are 
free from classical and quantum instabilities [91 do]. One stringent constraint on these models 
is imposed from the avoidance of hfth force carried by the extra scalar degree of freedom 
[m dH da], and a few successful models have been constructed m da [T6] . These models 
produce observationally consistent accelerating expansion preceded by the matter domination 

ra- 

Lately, the curvature singularity occurring at cosmological time scales have been found to be 
a serious problem in f{R) models [18]. It is well known that f{R) theories can be considered as 
equivalent scalar-tensor theories. It is easy to visualise curvature singularity in a f{R) theory by 
looking at the form of the potential appearing in its equivalent scalar-tensor theory in the Jordan 
frame [I9]. Due to the nonlinear motion of the scalar held, the oscillations around the potential 
minimum can make the held displaced to the singular point. The presence of the matter makes 
the occurrence of singularity more probable [18] . The hnite-time singularity in modihed gravity 
is described in [2ni[2T], whereas the singular behaviour of curvature in a contracting universe has 
been analysed in [22]. It is shown that past singularities may be prevented for a certain range 
of parameters. These singularities may also occur in future and can be avoided for hne-tuned 
initial conditions [iniEa]. It is also realised that the curvature singularities can be eliminated 
by adding an extra curvature term to the Lagrangian [22] |24]. The curvature singularity can 
also be seen in an astrophysical object. In this case, the singularity is analysed for suitable 
f{R) models applied to dense objects undergoing contraction in the presence of linearly time- 
dependent mass density [25l[26]. It is seen that the singularity is reached in a time that is much 
shorter than cosmological time scale. A detailed study of this issue can be found in [27] . 
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In this work, we study the issue of curvature singularity in a f{R) model proposed in [28] . 
We give special emphasis on the role of external matter that makes the effective potential 
shallower compared to the pure vacuum case. We find that it is impossible to find a parameter 
space where both the hfth-force constraint and the curvature singularity issues are resolved. 
This paper is organised as follows: In Sec. |2l we review the curvature singularity and fifth force 
constraints for a general function f{R)- The occurrence of curvature singularities and fifth 
force constraint in a specihc model of [28] is analysed in Sec [3l We extend the same analysis 
to a specific limit of the above mentioned model in Sec. 01 Finally, we summarise our results 
in Sec. [5l 


2 Curvature Singularity and Fifth-Force Constraints in 
f{R) Gravity 


In f{R) theories of gravity, the Einstein-Hilbert action is modified by replacing the Ricci scalar 
with an arbitrary function of R, 

S = T / d*x^f(R) + Sm (1) 

with tR = SvrG = Mpi, where Sm contains the matter degrees of freedom which does not include 
contributions from dark energy. For our convenience, we will write the arbitrary function in 
the following form f{R) = R + R^R), so that F{R) captures the modihcations of Einstein 
gravity. The effect of F{R) in the cosmological dynamics is relevant when local curvature 
becomes smaller than a characteristic infrared modihcation scale i?*. Above this scale, the 
gravity behaves approximately as Einstein gravity. 

By varying the action with respect to the metric tensor one can obtain the held equation, 
and taking trace of that held equation we arrive at 

3nFH(^) -2F-R + RF^niR) = k^T , (2) 


where T is the trace of the stress-energy tensor coming from the matter Lagrangian, and a 
comma with the subscript corresponds to the partial derivative with that quantity. The scalar- 
tensor representation of f{R) theory can be a found by identifying the term F^n as a dynamical 
degree of freedom, and in literature, it is dubbed as the ‘scalaron’ held 0 = Fr = — 1 [10] . 

The scalar curvature of the space-time is governed by the local value of the held (p. 

The above trace Eq. (I2|) can be conveniently written as 


dVj 


where 


dp 3 ^^ 


2F-RF 


,R) 


(3) 


(4) 
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It may not be always possible to invert the relation 0 = and obtain R{4>). Therefore it is 
convenient to write Vj in a parametric form 


dR 


dVj d(f) 
d(j) dR 


-{R + 2F — RF^r)F^iir . 


(5) 


Here, the subscript ' J' denotes the potential in the Jordan frame. The dynamics of the held 0 is 
governed by its potential Vj originating from the modihed gravity, and a force term that is pro¬ 
portional to the stress-energy tensor T. Integrating this expression we can hnd Vj(i?(0(x, t))). 
Thus, we can effectively describe the dynamics of modihed gravity by the dynamics of a scalar 
held whose value is uniquely determined by f\R). 


2.1 Curvature Singularity: 


The potential Vj{(j)) will have a global minimum at (pmin where cosmological evolution happens. 
But, as we will see, there is also a point in the held space, denoted by (psing, where the curvature 
scalar R diverges to inhnity resulting into curvature singularity. Typically, the points (pmin and 
4’sing would be Separated by a hnite held value with a hnite energy barrier. Therefore, the held 
can potentially reach to the singular point in the process of having small oscillations around 
its minimum, in particular when we consider objects with growing mass density [18]. The 
minimum of the potential (pmin can be obtained by equating dVj^cp) /dcp = R + 2F — RFr = 0 
which corresponds to a constant curvature solution DF^ji = 0 for the vacuum. Thus, (pmin is 
also a de Sitter point; see Eq.([2|). 

The appearance of the singularity can be seen both in the cosmological background [18] , as 
well as in the astrophysical dense object going under spherical collapse [25]. The Eq.([2]) can 
also simply be written as 


□0 


dVf 

d(j) 


( 6 ) 


where ^ J- (k^/3)T, i.e. Vj^ incorporates the ehects of matter. Now, we would like to 

study the dynamics of the scalar held in an astrophysical system whose mass density increases 
with time, i.e a collapsing object with time dependent With the assumption of weak 

gravity, the covariant derivative can be replaced by the usual hat space ones [27] . Again, in the 
approximation of isotropic matter distribution (spherical collapse), spatial derivatives can be 
neglected [25]. Thus, if we consider objects whose mass density changes with time T = T{t) in 
a homogenous way, the above equation simplihes to 


^ dVf{(p, t) ^ 

dt^ d(j) ^ ’ 

This is an oscillator equation with a time dependent nonlinear potential. Even if the initial os¬ 
cillation amplitude is small, due to the nonlinear behaviour of the motion, the above mentioned 
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singular point can be reached during the evolution of the held. As we will see, in many viable 
f{R) models when a collapsing astrophysical object is considered, 0 evolves to the singular 
point in a hnite time which is much shorter than the cosmological time scale ft; ~ 4 x 10^^ sec. 
This makes the theory not viable. 


2.2 Fifth Force Constraint: 


It is clear from the above discussion that the f{R) gravity theory contains a scalar degree of 
freedom in addition to the usual graviton. The effect of the scalar degree of freedom in the 
matter sector is evident when the theory is rewritten in the Einstein frame, where the gravity 
part is Einstein Hilbert type. Taking the conformal transformation on the metric 

= f,R , ( 8 ) 


the action in ([T]) transforms to 


S = 



1(VV.)" - VeW + 


(9) 


where, Ve is the scalar held potential in the Einstein frame 




Rffi - f 


( 10 ) 


All the quantities having tilde are dehned in Einstein frame. The mass of the scalar held ip 
can be as light as the present Hubble constant, and as a result, there will appear a long range 
hfth force mediated by the held. In order to satisfy the local gravity constraints one must have 
some screening mechanism which can screen the hfth force on the surface of the Sun/Earth. 
One such mechanism is known as chameleon mechanism na. The scalar held will be heavier 
in a dense object like earth, and light when the background density is less like the interstellar 
space. As long as the scalar held is heavy inside the surface of the earth, it will be frozen at 
the minimum of its ehective potential 


= Ve(V') + e“>^p , 


( 11 ) 


where p is the energy density of the matter held and cannot contribute to the outside held. 
The only contribution to the outside held can be from a very thin shell around the surface of 
the earth. The thin shell parameter is given by [12] 

P^out /' 10 ^ 

^ “ V&^test ’ ^ ^ 

where ipin, tpout are held values corresponding to the minimum of the ehective potential inside 
and outside of a test body (Sun/Earth), ^test is the gravitational potential on the surface at 
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the radius Tc of the test body and it is given by ^test = Mtest/^iifc where, Mtest = {4:7i/3)rcPin- 
Here pin is the average density of the test body. Arc is the width of the thin shell around the 
surface of the test body. In the usual approximation of -C 1 and F R, it can be shown 
that 'ipin ipout for all practical purposes, and the thin-shell condition reduces to [n] 

< VG^test—r^ (13) 

Tc 


^ j 5.97 X 10 (Solar system test), 

^ I 3.43 X 10“^^ (Equivalence Principle test). 

Therefore, in practice, for a particular model of f{R) gravity, we just need to calculate xfjout to 
see whether it satishes the fifth-force constraints mi, n. In the next section, we will discuss 
the curvature singularity problem for a particular form of f{R) gravity model in combination 
with the fifth-force constraint. 


3 A General Model 

In this section, we will consider a model proposed in [28] : 

f{R) = R + aR,l3\[l + {R/R,T ]~^'^, (15) 

where n, /9, a and i?* are positive parameters of the model, i?* is taken to be of the order 
of present day average curvature of the universe. The above function satishes the condition 
f{R) — )■ 0 as i? —0 and f{R) R + constant, for large R. This is necessary for the 
correct GR limit in early cosmological epoch m. Viable matter era demands n > 0,/9 > 0 
or u > 0, /S < —n mi. In this parametrisation, Starobinsky model corresponds to n = 2 [IT] , 
whereas Hu-Sawicki model can be obtained by plugging (3 = 1 [I5]. For the case of n > 0 and 
/3 < 0, the model becomes indistinguishable from the ACDM model [13]. In the original paper 
of [28], the model was analysed in the limit (3 ^ oo and for n = 1, and it was shown that the 
model is free from curvature singularity issue. We will discuss about the model in this limit 
in the next section. In this section, we will analyse the model in its full generality given by 
Eq. ([ISD. 

First, we briehy discuss the hfth-force constraints before moving into the curvature singu¬ 
larity issue which is the main focus of the article. In the limit of large curvature i? ^ i?*, the 
form of f{R) is reduced to 

/(R) = R + aR*/3|(R/R*)""/^ - l} . (16) 

As discussed in the previous section, to evade the fifth-force, the canonical scalar field in the 
Einstein frame outside the test body must satisfy the constraint of Eq. flTTl) . In the large 
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i/'out < 3.43 X 10 for ;S=1,5 and 10 
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Figure 1: parameter space of n and for different values of (3. Pink region shows allowed 

region for (3 = 1, purple region for (3 = 5 and yellow region for /3 = 10. 



curvature limit, the scalar held at its minimum of the potential outside the test body is given 
by [23] 

\/6 


= 1.1(1+ FR)~^r’.« = 


2k 


-na 


9 \ ---1 

^ pout ' ^ 

R* 


(17) 


where we have used the minimisation condition of Ves in evaluating R = K?Pout at the minimum 
of the effective potential. If Ri is the curvature at the de Sitter minimum Ri = K^Pcrit, we 
dehne a dimensionless variable Xi = Ri/R^ to obtain 


\'^ out 


Vq 

~ —na 
2k 


f ^Ipout 
V Pcrit 


---1 


(18) 


Now pcriticai — 10“^®gm/cm^, and pout — lll~'^‘^gm/cm? (typical baryonic/dark matter density). 
Because of F{R) <C f{R), ct should be of order of unity. In this case, the hfth-force constraint 
of Eq. flT^ can be easily translated to the condition n/(3 >2 for xi ~ 0{1). For general values 
of parameters the allowed regions are shown in Fig. [H where we have used the condition of 
de Sitter minimum in writing a in terms of other parameters. It is clear that the hfth force 
constraint is satished when n/(3 > 2 [23]. It can be been seen that in the limit of (3 ^ oo, 
the model immediately violates the hfth-force constraints [23], even though it can evade the 
curvature singularity problem [28] . 
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Figure 2: Vj/R^, vs. 0. In both figures a = 2.0. In the left figure, n = 2 and (3 = 1,5 ,10 and in 
the right figure n = 1, 2, 3 and (3 = 1. The minimum of a potential is marked with a dot sign. 

3.1 Curvature Singularity: Static Analysis 

In this section, we will discuss the issue of curvature singularity in the Jordan frame for the 
model described by Eq. (1T5|) . Looking at the form of the potential Vj, we argue for the presence 
of curvature singularity in this model. A preliminary analysis was done in [23], but in the 
Einstein frame. Moreover, the true nature of the singularity can be understood when the 
effects of the matter is taken care properly via [18]. Typically, the existence of matter 
pushes the minimum point close to the singular point. We will show it for ‘Log model’ in the 
next section - See Fig. HI In the next subsection, we will analyse the dynamics of the field by 
solving the equation of motion to show that the field indeed reaches to the singular point in a 
time scale much smaller than the age of the Universe. 

For our case, the scalar degree of freedom 0 in the Jordan frame can be identified as 

cl> = -na{R/R,T-^[l + {R/R,T]-^'^-^ . (19) 

Note that R ^ oo when 0 —)■ 0, i.e fsing = 0 in this model. At the same time, we also note 
that positive values to 0 correspond to negative curvature except for n = 1 case. In Fig. E] 
we plot the Jordon frame potential Vj 0. The potential is multivalued which is common to 
many f{R) gravity models, and in this case the physically relevant potential is the lower one. 
As described in the previous section, the de Sitter minimum of the potential, about which 
scalar field oscillates, can be obtained from the condition dVj/dcj) = 0 which is equivalent to 
the condition dVE/dip = 0. We see from the plot that as we increase the value of (3, the height 
of the potential barrier between the de Sitter point and the singularity point increases, and 
also the de Sitter minimum of the potential shifts away from the singularity in the field space. 

^The analytical expression is complicated, and is not very illuminating. Thus, we do not show it explicitly 
here. 
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Thus, for a given n, larger values of /9, the de Sitter minimum is at safer distance from the 
curvature singularity. This is in accordance with the hndings of |23] (in Einstein frame) where 
in the limit of /5 —?■ oo, the potential barrier becomes inhnitely large, thus inaccessible by the 
field and effectively hiding the singularity [2^ 128] . On the other hand, looking at the right 
panel plot of Fig. |2]we see that making n large pushes the minimum close to the singular point. 


3.2 Curvature Singularity: Dynamical Analysis 


We consider a self gravitating system (e.g astronomical objects of dark matter cloud) whose 
mass density is changing with time. This is in contrast to the static conhgurations in f{R) 
theories [291 EDI EH. We solve Eq. (IT)) for the contracting homogeneous and isotropic cloud of 
pressureless dust whose density pm is much greater than the critical density pcrit- For example. 
Pm ~ 10“^^ gm/cm^ for a dust cloud in a galaxy. We parametrise the trace of energy-momentum 
tensor T as 

T{t) = —Tq{ 1 + t/tch) , ( 20 ) 


where, Tq = pm is the energy-momentum density at t = 0, and tch is the characteristic time for 
density variations. It can be estimated by tch ~ d/v, where d would be the typical dimension 
of the collapsing object and v is its velocity during collapse. Here, we assume that the object 
collpases only under the effects of gravity and therefore velocity v is nothing but the escape 
velocity on the surface of the object. Consequently, tch can be calculated by tch ^J?>/{%'KGpm) 
which comes out to be ~ 1.34 x 10^^ sec for pm ~ cm^. The time within which the 

system meets curvature singularity is denoted by tsing- For a physical system, tsing must be 
smaller than the age of the Universe denoted by ~ 4 x 10^^ sec. We have assumed that the 
gravitational held of the object is weak and hence the background metric is a hat Minkowski 
metric, and because of the homogeneity and isotropy D’Alembertian has been replaced with 
ordinary time derivative [27] . 

With the above mentioned variation of matter density, the governing equation for the scalar 
held looks like 

-|- 2F + R — RF ji — k^Tq (1 -|- t/tch) = 0 . (21) 

Changing the variables from R and t to y = p{R^/R) and r = t/tch we obtain 


/2 

g ny _ 2 

0 y 


.-n/y 


(1 + r) 


1 2«/3 

y~^ + ^ 

V 


+ nap ^ ^ 


n J 


( 22 ) 


where, p = and Tch is given in terms of model parameters and characteristic time tch 


Tch = { [k^To/R^) 


2 + 


^ (9 R* 

n-\- (5 3na 


ich • 


(23) 


In deriving the above equation of motion, we have used the Eq. (ITbU . The singularity is reached 
when y = k^Tq/R becomes zero. The nonlinear behaviour of Eq. (|22D can amplify the initial 









Figure 3: The oscillation amplitudes are plotted against physical time for several parameter 
choices of the model. k^Tq/R vs. t. In both the figures, tch = 1-34 x lO^^sec, 77 = 10^ and 
a = 2. 


n 

13 

nIP 

(Sgc) 

2 

1 

2 

8.6 X 10^° 

2 

2 

1 

6.9 X 10^5 

4 

5 

0.8 

2.5 X 10^^ 

2 

5 

0.4 

1.0 X 10^5 

2 

100 

0.02 

3.9 X 10^5 

2 

10^ 

2 X 10-5 

4.5 X 10^5 


Table 1: tgi^g for different parameter values of the model dehned by Eq. ffTSj) . We have chosen 
a = 2. For all these choices of parameters, the time scale of singularity is much smaller than 
the age of the Universe. 

oscillation amplitude around the potential minimum and the held can reach to the singular 
point in a time scale that is smaller than the age of the Universe. The equation has been solved 
for 7/'(0) = 0 and |/(0) = 1 initial values, and it is consistent with R/R^, ^ 1 approximation. 
The solution of above equation for tch = 1-34 x lO^^Fec and different values of n and (3 is given 

in Fig. [3l 

It can be seen that the amplitude of the oscillations increase with time, and after a hnite 
time the value of y approaches zero which is a singular point. It is necessary that the time 
scale within which y goes to the singular point must be shorter than tu. As an example, taking 
i?* ~ 1/t^, tsing = 2.5 X 10^“^ Sec for n = 4, = 5 and a = 2. It is clear from the plots 

that higher values of (3 and lower n values are needed for the considered model to remain free 
from curvature singularity. This is consistent with the primary idea that we had obtained by 
looking at the potential Vj(0) for the scalar held. The time tsing has been listed in Table. [Ufor 


9 

























different values of parameters n and jS. It is evident that the relevant parameter for the time 
scale of singularity tging is the ratio between n and /3. The singularity time tsing increases with 
smaller values of n/(3. 

Thus, it can be concluded that singularity can be avoided in principle by taking very large 
values of (3 [28] and/or small values of n. On the other hand, as discussed earlier, local gravity 
tests constrain the values of the parameters n//3 > 2. Therefore, it is nearly impossible prac¬ 
tically to hnd a region in the parameter space where the model of Eq. flT5|) can be free from 
curvature singularity, as well as it can evade the local gravity tests. 

4 ‘Log’ Model 

In the last section, we have seen that the curvature singularity can be avoided in principle by 
taking large f3 and small n limit of the general model given by (1T5|) . In [28], it has been claimed 
that if we take n = 1 and f3 ^ oo limit of Eq. ffTSj) . the generalised model is reduced to the 
following logarithmic function: 


f{R) = R — aR^, In (1 -I- R/R^) . 


(24) 


Subsequently, this particular form of the model has been analysed and it is shown that the 
curvature singularity is absent [28], even though it violates the hfth force constraint of n//3 > 2. 
First of all, this reduction is only possible under the assumption i? -C i?* and it breaks down 
immediately for the analysis of curvature singularity. Therefore, we consider Eq. (j24|) as an 
independent model and revisit the issue of curvature singularity in the presence of matter. 

The scalar held (j) in the Jordan frame for the model can be expressed as 


a 


1 J- R/R^f 


(26) 


and we see that 0 —)■ 0 corresponds to inhnite R corresponding to the curvature singularity. To 
investigate the issue, we hrst hnd out the potential 


Vj{(j))/R^ = -^0- + ~ + ^aln(-l/(/)) , (26) 


and note that Vj{(j)) —?■ cxd as 0 reaches the singularity point. Naively, the presence of the inhnite 
potential barrier was the argument in claiming the absence of curvature singularity [23], [28] . 
But including the matter source in term of a collapsing cloud of Eq. ([20]), the ehective potential 
can be written as 



1 1^2 2 ^ 1 , 




For diherent values of k'^Tq/R^ at t = tch, the plot of vs. 0 has been shown in Fig. [4]. 
The de Sitter points are marked by stars in the plot. From the hgure, it can be seen that 


10 







Figure 4: Left: Vj® vs. 0 where dots represent the potential minimum. In the hgure, k^Tq/R^ = 
0,10,50 and a = 1.5. Right: k^Tq/R vs. t/tch for rj = 10^, 10^, 10® with a = 1.5. 


the minimum of the potential shifts towards the singularity (0 = 0) as the density of matter 
increases [H] . One can conclude that the possibility that the scalar held (j) may hit the curvature 
singularity is more likely in the presence of matter than in the vacuum. But, it requires the 
understanding of dynamics of the held which we will analyse now. 

We investigate the presence of curvature singularity by solving Eq. 0 numerically in con¬ 
tracting astrophysical objects as we have done in the previous section. The equation of motion 
in this case reads 


y 


T, 


ch 


2ln{y) — 2ln{ri) -|- —- -|- 1 — — (1 -|- r) 
ay a 


= 0 , 


(28) 


where the approximations at large curvature limit F{R) — —aRJn {R/R^) and ~ —{aR^)/R 
have been used. Here, Tch is given by 


(29) 

The solution of the above equation |/(r) is plotted in Fig.lHfor pm = 10“^^, 10“^^ and gm / cm^ 

i.e. for T] = 10^, 10^ and 10®. The characteristic time tch for these densities are 1.34 x 10^^, 1.34 x 
10^^ and 1.34 x 10^^ sec correspondingly. The initial conditions have been chosen as ?/(0) = 1 
and |/'(0) = 0. It is seen that the oscillations of y grows with time and eventually meets the 
singularity. It is evident from the hgure that the singularity is reached earlier for larger r]. 

It was noted earlier that the model is singularity free in the absence of matter because of its 
large potential barrier at the singular point. But we show that in the presence of astrophysical 
densities in a collapsing object the system indeed can reach to the singular point. 
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5 Discussion and Conclusion 


In this work, we have investigated the existence of curvature singularity problem in f{R) model 
given by Eq. flTSl) and by Eq. fl2^ |28j . For our analysis we have chosen this particular model as 
this form of f{R) incorporates several f{R) models well known in the literature. We have made 
the dynamical analysis of the fluctuations of the associated scalar field around the minimum of 
the potential to show that the held indeed reaches to the singular point when the effects of the 
matter are considered. 

We have carried out our analysis of curvature singularity in the Jordan frame where the 
theory has been rewritten in terms of a scalar held (p. The singularity is reached when 0 = psingi 
and in our case, psing = 0. The scalar held oscillates around its de Sitter minimum of the 
potential Vj. The potential barrier at the singularity point, and the distance between the de 
Sitter point and the singular point are hnite. Therefore the scalar held can reach to (psing 
during its cosmological evolution or during the collapse of dense astrophysical objects. In this 
work, we have shown this ehect by solving the nonlinear evolution equation of the held. With 
increasing value oi n/f3 in the considered model, the system reaches to the singularity earlier. 
In particular, we have solved the stress equation for collapsing astrophysical object and shown 
that the Ricci scalar R oscillates and after a time t ~ tging, it reaches to a divergent value. 
The time within which singularity is reached, i.e. tging, is less than the age of the Universe for 
generic values of n and /J. On the other hand, local gravity tests put the constraint on the 
model parameters such that n/(3 > 2. Taking large enough values of 0 can in principle help us 
avoiding the singularity in cosmological time scale, but it immediately violates the hfth force 
constraints. Therefore, we conclude that the model can not satisfy constraints coming from the 
curvature singularity and hfth force simultaneously. 

On the other hand, the ‘Log Model’ as discussed in Sec. 0] seems to be free from curvature 
singularity [23],[28]. But, this conclusion was drawn in the absence of matter. When we 
incorporate the effects of matter in the calculations, we see that the minimum of the effective 
potential moves closer to the singular point and the height of the barrier at psing becomes 
smaller. In fact, the solutions of the scalar held dynamics show that the singularity is reached 
in a time scale smaller than the age of the Universe. 

Looking at the present analysis in conjunction with previous understanding, it is tempting 
to think that it is very difficult, if not impossible, to construct any function f{R) that simul¬ 
taneously satishes both hfth force constraint and resolves the singularity problem. As we have 
seen, and also has been noted earlier [23], the height of the potential barrier goes opposite to 
satisfying the hfth force constraint. For example, in the considered model, the height of the bar¬ 
rier (thus avoidance of singularity) is proportional to 0, but having large (3 immediately violates 
the hfth force constraints. Therefore, it would be very interesting in hnding generic structure 
of the function f{R) that can simultaneously satisfy both the above mentioned constraints. We 
hope to look into this issue in future. 
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